Stability of the U(l) spin liquid with spinon Fermi surface in 2+1 dimensions 
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We study the stability of the 2+1 dimensional U(l) spin liquid state against proliferation of 
instantons in the presence of a spinon Fermi surface. By mapping the spinon Fermi surface into 
an infinite set of 1+1 dimensional chiral fermions, it is argued that an instanton has an infinite 
scaling dimension for any nonzero number of spinon flavors. Therefore, the spin liquid phase is 
stable against instantons and the non-compact U(l) gauge theory is a good low energy description. 



I. INTRODUCTION 

Fractionalizcd phase is a novel state of correlated 
many-body systems where low energy excitations carry 
fractional quantum numbers of microscopic degrees of 
freedom. In 1+1D, a spin-charge separation, which is an 
example of fractionalization, can naturally occur due to 
the low dimensionality [l|. In 2+1D, fractional quantum 
Hall states support excitations which have fractional elec- 
tric charges[2j. Finding a fractionalizcd phase in time- 
reversal symmetric 2+1D systems is an outstanding prob- 
lem in condensed matter physics 0, 3- 

In fractionalizcd phases, there exist non-local correla- 
tions which are not captured by the conventional sym- 
metry breaking picture [5(. Those correlations are asso- 
ciated with a condensation of stringy objects in space[6| 
or membranes in space-time|7(. It turned out that the 
most natural framework to describe those correlations is 
gauge theory, where the gauge field describes transverse 
fluctuations of condensed strings or membranes. 

Spin liquid is a fractionalized state where an elemen- 
tary excitation is spinon which carries spin 1/2 but 
no chargej3j. Among a variety of possible spin liq- 
uid states [5j, the state which has fermionic spinons and 
an emergent U(l) gauge field has been proposed for 
many 2+1D strongly correlated electron systems includ- 
ing high temperature superconductors, frustrated mag- 
nets and heavy fermion systems. Although high T c 
superconductors have the conventional superconduct- 
ing ground state, the normal state shows non-Fermi 
liquid behaviors which are possibly due to a proxim- 
ity to a spin liquid state [J, [IJ]. Frustrated magnets 
are simpler systems than the high T c cuprates in that 
there is no low energy charge mode. At the mo- 
ment, there are promising candidate materials |l2j, ll3| 
for which the U(l) spin liquid states with fermionic 
spinons have been proposed[8|, [9|, |lOj. Related fraction- 
alized phases have been studied in heavy fermion sys- 
tems near magnetic quantum critical points [lj] and frus- 
trated bose systems J15J. In the U(l) spin liquid states, 
fermionic spinons have either nodal points or Fermi sur- 
faces. The gapless spinons are strongly coupled with the 
U(l) gauge field at low energies and there is no well- 
defined quasiparticleQJ, G3, Qj, 01, M, S3 . 

Because of underlying lattice structures, the U(l) 
gauge field is compact, which allows for a topological 



defect called instanton (or monopole). Instanton, as a 
localized object in space-time, describes an event where 
the flux of the gauge field changes by 2ir. Understanding 
dynamics of instantons is crucial because the fractional- 
ized state is stable only if instantons are suppressed in 
the long distance limit. It has been known that if there is 
no gapless spinon, instantons always proliferate, result- 
ing in confinement. In this case, spin liquid states are not 
stable and spinons are permanently confined[22|. In the 
presence of gapless spinons, it is possible that the gauge 
field is screened and instanton becomes irrelevant in the 
low energy limit. If this happens, fractionalized phase is 
stable and spinons arise as low energy excitations. 

If there are a large number of gapless spinons which 
have the relativistic dispersion near nodal points, it has 
been shown that instanton is irrelevant at low energies 
and the fractionalized phase is stable[23|]. However, it is 
largely unknown whether the spin liquid phase is stable 
in the physical cases where the number of spinon fla- 
vors is relatively small. In the presence of spinon Fermi 
surface, it has been speculated that the abundance of 
low energy spinon modes may stabilize the fractional- 
ized phase more easily. However, dynamics of instan- 
tons in the presence of non-relativistic spinons has not 
been well understood. There have been several RPA 
studies[H, [H, [H, H3, IH IH, but currently there exists 
no non-pcrturbative analysis on the fate of instantons in 
the presence of spinon Fermi surface. Particularly, a lack 
of the conformal symmetry makes it hard to treat the 
problem in a non-perturbative way which is required be- 
cause instanton itself is a non-perturbative phenomenon. 

In this paper, we provide a non-pcrturbativc argument 
which supports the idea that the U(l) spin liquid state 
with spinon Fermi surface is indeed stable against pro- 
liferation of instantons for any nonzero N, where N is 
the number of spinon flavors. The paper consists of the 
three parts. In the first part (Sec. Ill), we ignore fluctu- 
ations of the non-compact component of the gauge field 
and calculate the scaling dimension of instanton at the 
fixed point described by free spinons. To do this, in Sec. 
II, we formulate low energy modes near the Fermi surface 
in terms of an infinite number of 1+1D chiral fermions. 
Since an instanton is a localized source of 2-7T flux in space- 
time, the fermions which move in 1+1 dimensional sub- 
spaces have to enclose the half of the solid angle around 
the instanton and acquire phase ix, when they are trans- 



ported around the instanton at a sufficiently large dis- 
tance. This is illustrated in Fig. [21 Therefore, an in- 
stanton operator corresponds to a twist operator of the 
1 +1D chiral fermions. The scaling dimension of an in- 
stanton is infinite because there are infinitely many 1+1D 
fermions parametrized by the direction of their velocities 
(or angular momentum), and each fermion contributes a 
finite scaling dimension to the total scaling dimension of 
the instanton operator. In the second part (Sec. IV), 
the fluctuations of the non-compact gauge field are con- 
sidered together with instantons. To control the gauge 
fluctuations, we consider a large N limit. In this case, 
vertex corrections are negligible and we can obtain a def- 
inite scaling transformation under which the low energy 
theory remains invariant. The key difference from the 
previous studies [19|,|2l| is that in the present approach all 
points on the Fermi surface are treated on the equal foot- 
ing rather than focusing on a local patch in the momen- 
tum space. This enables us to define the scaling dimen- 
sion of the instanton operator, taking into account the 
whole Fermi surface. With the fluctuating non-compact 
gauge field, fermion modes which have different Fermi ve- 
locities arc no longer decoupled, and we can not simply 
sum the scaling dimensions of different modes as we did 
in the non-interacting case. However, in the low energy 
limit, only small angle scatterings are important because 
momenta of the gauge field are scaled down while the 
circumference of the Fermi surface is unchanged under 
the scale transformation. This implies that two fermion 
fields on different points on the Fermi surface arc essen- 
tially decoupled at low energies. Therefore, there are 
still infinitely many independent 1+1D fermion modes 
which contribute to the scaling dimension of instanton 
at low energies. By using this property, we can argue 
that the scaling dimension of an instanton is infinite at 
the interacting fixed point too. Finally, in Sec. V, we 



consider the case with a small N of the order of 1 which 
is directly pertinent to the U(l) spin liquid state with 
two flavors (spin up and down) of spinous [a [j| proposed 
for k - (BEDT - TTF) 2 Cu 2 {CN) 3 \v^. With a small 
N, the Fermi surface is strongly coupled with the fluc- 
tuating gauge field and vertex corrections can not be ig- 
nored. This makes it difficult to find an explicit form of a 
scaling transformation for the strongly interacting fixed 
point. However, one can see that the essential properties 
which make the scaling dimension of instanton infinite 
does not depend on the specific form of a scaling trans- 
formation. Actually, the existence of an extended Fermi 
surface and the fact that only small angle scatterings are 
important at low energies are enough to argue that the 
scaling dimension of instanton remains to be infinite and 
instantons are irrelevant at the strongly interacting fixed 
point for any nonzero N. 

II. ANGULAR REPRESENTATION OF FERMI 
SURFACE 



We start by considering N flavors of fermions coupled 
with a compact U(l) gauge field in 2+1D, 
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Here 'F,- is the fermion field with N flavors, j = 1, 2, .., N 
and a M = (ao, a) is the U(l) gauge field with \x = 0, 1, 2. 
Up is the chemical potential and g, the gauge cou- 
pling. ffj, v is the field strength tensor. Summation over 
the repeated flavor index j is implied. In the energy- 
momentum space, the action becomes 
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Here p, I denote energy-momentum vectors and e p = -L- 1 [ip . Integrating out high energy fermion modes outside a 

momentum shell with a width A near the Fermi surface, we obtain the low energy effective action S = So + Si, where 

S = J diodkdd (iuj + k) if;* (u,k, 6)i>j (uj 7 k,6) 
1 
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Here, the Fermi velocity has been set to 1. Fermion mo- 
mentum is represented in the polar coordinate [3fJ where 
k = |k| — kp is the deviation of momentum from the 
Fermi surface in the radial direction and 9 is the an- 
gular coordinate as is shown in Fig. Q] (a). We use 
the approximation, J dk = J d\k\ |k| J dd sw kp J dk J d9 
and redefine the fermion field as ipj(u>,k,6) = 
kp /2 ^j (u>, fci = (k F + k) cos 6, k 2 = (k F + k) sin 8). K ~ 
Nkp is the diamagnetic term, ag — kg ■ a is the spa- 
tial gauge field parallel to the fermion momentum along 
kg = (cos 0, sin 6). qi = kg ■ q and q t = (kg x q) z are 
the momentum components of the gauge field which are 
parallel and perpendicular to kg respectively. Note that 
a o(^> lii 1u 6) an d oe( I/ , Ih Qt\ &) m t ne second line of Eq. 
© implicitly depend on 9 because qi and q t are measured 
with reference to kg as is shown in Fig. Q] (b). A is the 
momentum cut-off of the fermions near the Fermi surface 
and A is the cut-off of the gauge field. For A « kp, we 
can ignore the quadratic term k 2 /2m which is irrelevant 
at low energies. 



mentum quantum number, 
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In the 1+1D real space, the action for the chiral fermions 
becomes 



S = ^2 drdx q Pjn( T ' X ) 



Or 



i(d x - ia 8 g{T,x,x t = n/kp)) ip jn (T,x), (5) 



where the 1+1D gauge field 'felt' by the chiral fermions 
is given by 



a s g{T,x,x t ) 
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III. FREE FERMIONS 

The gauge field can be decomposed into the singu- 
lar part which includes instanton configurations and the 
non-singular part which describes fluctuations of the non- 
compact component. First, we ignore the non-compact 
gauge field and examine the effect of instantons on the 
free fermions. For this, we consider a background gauge 
field a s generated by an instanton located at r = and 
x = in space and time. For the singular part of the 
gauge field, we use the temporal gauge where ag = 0. 

If the Fermi surface has the rotational symmetry, the 
field strength for single instanton centered at the ori- 
gin has the rotational symmetry too. This enables us 
to choose a 8 e (v, qi,qt] 0) to be independent of 9. Because 
of the rotational symmetry, it is convenient to introduce 
1+1D chiral Fermi fields which have good angular mo- 



Here x and xt represent the displacements from the in- 
stanton in the directions parallel and perpendicular to 
kg for some 9 respectively. Because of the rotational 
symmetry, ag(r, x, xt) is independent of 9 and the choice 
of 9 does not matter. The action in Eq. ((5]) describes 
an infinite set of 1+1D chiral fermions coupled to the 
background gauge field, o|(t, x, x% = n/kp). The chi- 
ral fermion with angular momentum n 'sees' the 2+ID 
gauge field projected onto a plane which is perpendicular 
to the xi — X2 plane and shifted by njkp away from the 
origin as is shown in Fig. [2] (a). The range of the mo- 
mentum integration in Eq. ([6]) is restricted to be within 
the strip with the width A as is shown in Fig. [T] (b) and 
ag(r, x,Xt) represents slowly varying configurations of the 
gauge field in space and time. Nevertheless, the compo- 
nents with large momenta become unimportant in the 
long distance limit, and a s g (T,x,Xt) accurately describes 
the true configuration of instanton far away from the cen- 
ter. For an instanton whose field strength is isotropic in 
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FIG. 1: Low energy modes in the momentum space (a) for 
the fermions and (b) for the gauge field. A and A are the cut- 
off momenta for the fermions and the gauge field respectively. 
In (a), the bold circle represents the Fermi surface and the 
arrow connecting the two filled circles represents a momen- 
tum transfer from the gauge field. The shaded strip in (b) 
represents the points of momenta which are included in the 
background gauge field felt by the 1+1D chiral fermions. 



space and time, the gauge field is given by 



a 2 



= 0, 



•('2 



2r(r — r) ' 

Xl 

2r(r — r) 



(7) 



with r = \Jt 2 + x\ + x\. In this gauge, there is a Dirac 
string stretched along the positive r axis. The presence 
of the Dirac string is not important because the infinitely 
thin tube of 27r flux can always be placed inside a halo 
of an underlying lattice and the unit flux can be gauged 
away. The gauge field in Eq. ([7]) represents an instanton 



FIG. 2: A 1+1D chiral fermion which has an angular mo- 
mentum n is coupled to a background gauge field projected 
to the plane which is represented as the vertical plane in (a). 
The arrowed circle in (a) is a trajectory of a fermion trans- 
ported around an instanton at distance L and the arrowed 
circle in (b) is the trajectory of the fermion projected on the 
unit sphere. The area of the shaded region in (b) is the solid 
angle enclosed by the fermion around the instanton which is 
represented as the black dot at the center of the sphere. In 
the large distance limit (L >> h/If), the fermion acquires 
phase 7T as it moves around the instanton. 



with the Lorentz symmetry. In the presence of the non- 
relativistic fermions, the space-time isotropy is lost and 
the field strength will be redistributed. With the broken 
Lorentz symmetry, the spatial rotational symmetry in the 
xi — X2 space may or may not be broken. In the following, 
we will first consider the case with the spatial rotational 
symmetry and then consider general cases without the 
symmetry. 

As an 1+1D fermion is transported within a plane at 
distance L from its origin as in Fig. d](a), it encloses the 



solid angle, 0, n (L) = 2tt I 1 — 



y/(Lk F ) 2 +n 

sphere around the instanton as is shown in Fig. [2] (b). 
Since an instanton is the source of flux 2tt, the fermion 
acquires a non-trivial phase $„(£). Without loss of gen- 
erality, we can defined the phase angle within the interval 
(— 7r,7r]. For the Lorentz symmetric instanton configura- 
tion in Eq. (J7J), the phase angle is the half of the solid 
angle and becomes 



,i ■, where <7j n (7r) is the operator which twists the boundary 



$ n (L) = sgn(n)ir 1 



y/(Lk F ) 2 +n 2 
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For non-isotropic configurations, <£> n (L) will be differ- 
ent from Eq. (jSJ). However, the explicit form of $ n (X) 
is not important for the following discussions. What is 
important is the fact that |$„(L)| — > 7r for any finite n 
as L — > oo in the presence of the spatial rotational sym- 
metry. This is because trajectories of fcrmions projected 
onto the unit sphere around the instanton will eventually 
follow a big circle for any angular momentum n in the 
large L limit, and the flux enclosed by any half sphere 
that cuts through the north and south poles is always 7r 
due to the spatial rotational symmetry. Therefore, in the 
long distance limit, an instanton twists boundary condi- 
tions of all fermions from the periodic condition to the 
anti-periodic one. Roughly speaking, at a length scale L, 
an instanton twists boundary conditions of the fermions 
which have angular momenta \n\ < kpL by it. 

Having understood that an instanton operator corre- 
sponds to a twist operator, we can determine the scaling 
dimension of instanton. To do this, we represent the 
1+1D space in terms of a complex variable, z = r — ix 
and rescale the fermion fields to write the action in the 
standard form, 



1 N oo p 

Sfree = — ^ ^Z / dzdI ^jn^Pj 



(9) 
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where d = ■§=. This free theory is invariant under the 
scale transformation, 

z = bz , 
z = bz , 
j, jn (br' ,bx') = b-^% n (r,x) (10) 

with 6 > 1. With instantons, the free theory is perturbed 

as 



S = S 



free 



dzdz 
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where a is the creation operator of an instanton or an 
anti-instanton and y is the fugacity. Since an instanton 
twists the boundary conditions of all fermions in the low 
energy limit, the instanton operator can be written as 



n, 



.W, 



(12) 



condition of tpj n by it. However, the scaling dimension of 
a is not well-defined because the operator ends up twist- 
ing the infinite number of fermions by the finite angle, w. 
Actually, not all fcrmions are twisted by the same angle 
at a finite length scale L 3l| . Therefore we introduce a 
regularized instanton operator, 



0"L = ILj, n <7jn(®n(L)), 



(13) 



where cr^ ($„(£)) is an operator which twists the bound- 
ary condition of tpj n by the angle $„(L). Physically, 
<jl creates the flux configuration near an instanton at 
a length scale L. Although a^ is not a true instanton 
operator, we can learn about the property of the true 
instanton by taking L — > oo limit of ctl- The point of 
introducing the regularized operator is that ctl is a well 
defined local operator which has a finite scaling dimen- 
sion, as will be shown below. 

Since the scaling dimension corresponds to the eigen- 
value of the scale transformation generated by the 
Noether current j v = x^T^, the scaling dimension d^ of 
the regularized instanton operator can be obtained from 



/dz 
— zT(z)a L (0), 



(14) 



where T(z) is the holomorphic energy momentum tensor. 
In the state-operator correspondence, we can view the 
scaling dimension as the 'energy' of the quantum state 
defined on the circle around the origin associated with 
the 'time' evolution in the radial direction. What the 
insertion of the ctl operator does is to twist the boundary 
condition of ij)j n by <&„ (L) . Then we can rewrite Eq. (fl~^ 
as 



dr = 



— *T{z) 



(15) 



twisted b.c. 



where we impose the twisted boundary condition for the 
fermion fields around the origin. Each fermion ipj n con- 
tribute a scaling dimension d($„(L)) = "p -J (see ap- 
pendix A for derivation) and the total scaling dimension 
for the regularized instanton operator becomes 



j— 1 n— — oo 



(16) 



Since <& n (L) approaches 7r as L increases, dj, diverges in 
the large L limit. For the isotropic instanton configura- 
tion, we have 



JV oo / , , 
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and it is easy to check that it diverges linearly with L. 
The renormalization group equation for the fugacity be- 
comes 



dyL 
din 6 



(2-d L )y L + 0(y 2 L ), 



(18) 



where y^ is the fugacity of the regularized instanton 
operator. The present approach does not allow us to 
calculate the higher order terms in y. However, from 
the linear term alone, we can readily see that for any 
N > the regularized instanton operator for a suffi- 
ciently large L (hence the true instanton operator de- 
fined as a = linii-jooiTi) is strongly irrelevant at the 
fixed point with y = 0. Namely, a small nonzero y will 
Sow to the fixed point with y = 0. Although the scal- 
ing dimension of the true instanton operator defined as 
d = lirn£,_ >00 dh is ill-defined (infinite) for any N > 0, the 
regularized scaling dimension diverges more rapidly with 
increasing L when N is larger. This is consistent with the 
physical intuition that the presence of more fermions re- 
sults in a larger scaling dimension of instanton via screen- 
ing. 

Until now, we have considered the case with the spatial 
rotational symmetry. If the symmetry is broken by an un- 
derlying lattice, the field configuration of an instanton is 
no longer symmetric under the spatial rotation. Here we 
will see that the conclusion reached for the rotationally 
symmetric Fermi surface holds in general cases too as far 
as the general Fermi surface can be obtained from the 
symmetric one through a smooth deformation. Perform- 
ing the Fourier transformations for the frequency and the 
radial momentum, we rewrite Si in Eq. ([3]) as 
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(19) 



For a general Fermi surface, ag has a non-trivial angular 
dependence on (9 + 9 )/2. To diagonalize the action, we 
have to use a different basis than the angular momentum 
basis. We consider the basis transformation, 



ipjn{T,x) = 



1 



dB f*(T,x,e)^(T,x,e), (20) 



(2tt) 1 /2 
where f n (T,x,9) satisfies the eigenvalue equation, 

f dO'a 9 (T, x, k F (9 - 9): (9 + e')/2)f n ( T , x, 9 ) 

(21) 



(2tt) 1 /2 

= a n (T,x)f n {T,X,6) 



and the normalization condition 

2~ I d0 fm( T , x >9)fn(T,X,0) = 5 n 



(22) 



One can always find such a basis because the kernel as 
satisfies the Hcrmitian condition, ag(r, x,k F (9 — 9 );(0 + 
6')/2) = a^{T,x,-k F (0-6');(9 + e')/2) at each r and x . 
In the new basis, the action for the 1+1D chiral fermions 
becomes diagonal, 



E 



drdx ipj n (T t x) 



Or 



-i(d x -ia n (T,x)) tpj n (T,x). (23) 

The phase acquired by the n-th fermion when the fermion 
is transported around the instanton is 



$„ = j> dxa n , 



and the total scaling dimension becomes 



(24) 
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where — n < $„ < n. For the rotationally symmetric 
Fermi surface, the index n represents the angular mo- 
mentum as before and $„ = 7r for all n in the low energy 
limit. As the Fermi surface (hence the field configura- 
tion of an instanton) is distorted, the distribution of $ n 
gets broadened around ±7r and d will decrease. How- 
ever, d can not change abruptly as the Fermi surface is 
smoothly deformed and d will remain infinite in the ther- 
modynamic limit unless there is a phase transition asso- 
ciated with the topology of the Fermi surface. Therefore, 
the scaling dimension of an instanton is infinite for a gen- 
eral Fermi surface as far as the Fermi surface is smoothly 
connected to the rotationally symmetric Fermi surface. 



IV. INTERACTING FERMIONS WITH A 
LARGE NUMBER OF FLAVORS 



Now we consider the whole theory by taking into ac- 
count fluctuations of the non-compact gauge field. For 
the non-compact gauge field, we will use the Coulomb 
gauge where V • a = 0[32j and drop the temporal gauge 
field which is screened out at long distances. In the fol- 
lowing, we will focus on the rotationally symmetric case. 
With the fluctuating non-compact gauge field, the inter- 
action term can be written as 




qsin 4> 



(b) 



(c) 



FIG. 3: Angular representations of (a) the fermions and (b) 
the gauge field. Note that the angle 8 for the fermions is 
measured from the X\ axis and the angle a for the gauge field 
is measured from the X2 axis so that the momentum of the 
gauge field becomes perpendicular to the Fermi momentum 
when a = 6. In (a), the arrow represents a process where 
a fermion changes its radial momentum and angle through 
an absorption of a momentum from the gauge field. In (b), 
cj> parametrizes a deviation of a from 0. In (c), the arrows 
represent the polarization vectors of the transverse gauge field 
at different points in the momentum space. 
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where a is the transverse gauge field. The momentum of the gauge field is also written in the polar coordinate as 
a(v, q, a) — a(is,qi = —qsin a, qi = qcosa). t a is the polarization vector of the gauge field with an angle a, where 
the angle is defined in such a way that t a becomes parallel to kg when a = 9 as is depicted in Fig. [3] We choose 
the polarization vector as t a = (cos a, sin a) for — tt/2 < a < ir/2 and t a = (cos(a + 7r),sirx(a + w)) for a > 7r/2 or 



a < — tt/2 as is shown in[3](c). In this definition, there arc discontinuities in t a at a = ±ir/2. However, this definition 
is convenient to make the reality of the gauge field explicit in the momentum space as a(~v, q,a + tt) = a*(y, q 1 a). 
Q(x) is the step function which ensures that the momenta of fermions lie within the shell of the width A near the 
Fermi surface. 

At low energies we have A ~ \fkph. because the Fermi surface is locally parabolic. Since kp >> A >> A, typical 
momenta of the gauge field are perpendicular to fcrmion momentum and much larger than A. Therefore the support 
of the <f> integration in Eq. (|2l))) is sharply centered at 4> = and (f> = tt with the width of the order of A/ A << 1. This 
allows us to use sin</> f=s ±</> and cosip sa ±1 in the low energy limit. Because both of the fermions with angles 9 and 
9 + tt are coupled with both of the gauge fields with angle 9 and 9 + tt, it is convenient to define two separate fields 
for the opposite-moving fermions and restrict the integrations of 9 to run from — 7r/2 to tt/2. At the same time, we 
allow q to run from —A to A and restrict the <f> integration to —A/\q\ < <j> < A/\q\ to implement the thcta functions 
in Eq. (|26p . Then the action can be written as S = So + Si , where 



/oo /"^V^ / 17r /^ 

duj dk ' d9 (iu> + sk) tp* js (lu, k, 9) ip js (u, k, 9) 
-oo J-A/2 J-tt/2 

/•oo /-A /-tt/2 / j \ 

+ / dv J dq\q\ J da f t^^ + Q 2 ) + K J a (^ Q, a)a(-v, -q, a), 

-i |.oo /-A/2 pjr/2 

5l = - E T^ZVi 



r A/2 

dv I dk I d9 

-A/2 J-tt/2 



fA r-A/\q\ 

dv I dq\q\ / 

-A J-A/\q\ 



(27) 



^ (2-) 3 / 2 
sa{u, q,9- <j>)il)* s (w + v/2, k + q<j>/2, 9 + sq/2k F ) ip ja (u - v/2, k - q<t>/2, 9 - sq/2k F ) , 



(28) 



r 



where s = 1,-1 labels fermion fields on the two sides of 
the Fermi surface at each angle —tt/2 < 9 < tt/2, defined 
as 

ip j+ (uj,k,9) = ij)j(u),k,9) 

V>j_(w,M) = ipj(u),-k,-sgn(e)Ti+9), (29) 

and a negative q of the gauge field represents the opposite 
momentum as 



a(v, —\q\,a) = a(v, \q\, a + tx) 



(30) 



with —tt/2 < a < tt/2. ipj + and ipj- have the opposite 
velocities and they form a two-component 1+1D Dirac 
fermion. 

The present approach is conceptually analogous to the 
bosonized descriptions of Fermi surface [33, [3J, [3a, [361 ] - 
where chiral bosons describe low energy particle-hole ex- 
citations near Fermi surface. However, there is an impor- 
tant difference. In the previous bosonized description of 
the Fermi surface coupled with the U(l) gauge field [3q. 
the Fermi surface is taken to be locally flat within each 
momentum patch, and the local curvature is not taken 
into account. The present formalism captures the local 
curvature effect of the Fermi surface, which is important 
to reproduce correct low energy behaviors [2l| . The key is 
to consider all points on the Fermi surface on the equal 
footing, not treating the Fermi surface as a sum of lo- 
cally flat Fermi segments. The way the curvature effect 
is implemented in Eq. (|2"8")l is explained in Fig. [U 

Here we assume that N >> 1 in which case the fluc- 
tuations of the gauge field are controlled. In the leading 



order of the 1/N expansion, vertex corrections can be 
ignored [1^. The fluctuating gauge field and the gapless 
fermions lead to singular self energies and the single par- 
ticle quantum effective action of the fermions and the 
transverse gauge field becomes (see appendix B) 



o — 



E 

s=±l 



dudkd9 



ic sgn(w)|w| 2/3 + sk 



■<Pjs(v,k, 9)tp js (ui,k,9) 



dvdq\q\da 



7 H 
7 M 



xq 



a*(v, q, oc)a(v, q,a), 

(31) 



where c, 7 and x are constants. Because of the singular 
quantum corrections, the scaling transformation in Eq. 
(TiT))) is no longer a symmetry. Instead, we have to rescale 
energy and momentum as 



U! 


= &-v, 


;/ 


= irV, 


A 


= b-^ 3 A, 


A 


= b~ l ' 3 K, 


k 


= b~ 2 / 3 k, 


<1 


= b- l 'U- 



(32) 



Note that the momentum of the fermion in the radial 
direction and the momentum of the gauge field should 
scale differently. If we apply this new scale transforma- 
tion to S\ in Eq. (|2"8")l . we readily notice that the action 




FIG. 4: The effect of the local curvature in the Fermi sur- 
face. As a fermion with a momentum A on the Fermi surface 
absorbs a momentum q which is tangential to the Fermi sur- 
face at the point, the fermion ends up with the momentum B 
which has a higher energy by q 2 /2kp due to the curvature of 
the Fermi surface. This can be seen from Eq. (|28|l as follows. 
In order for q to be tangential to the Fermi surface at the 
initial momentum, (j> in Eq. ([280 has to be q/2kF- Then the 
energy of the final state is larger than that of the initial state 
by q4> = q 2 /2kF- 



can not be made invariant unless the angular variables 
9 and <j) are rescaled as well. This is because momenta 
of the fermions and the gauge field mix with the angular 
variables through k + q<f> and 9 + sq/kp. To make the 
action invariant, we should assign the scaling dimension 
1/3 to the angular variables. This is an anomalous scal- 
ing dimension of the angular variables which arises solely 
from quantum effects. The whole action is invariant if 
we rescale 

9 = b- 1/3 e', 
M^^'.b-^k'x 1 '^') = b 4 ^^' a (uj',k,e'), 

a(&-V,6- 2/ V,6- 1/: y) = b i/3 a'(v',q\cj)') (33) 



to (— bir/2, bir/2). In the low energy limit where b — > oo, 
9 becomes a non-compact variable which runs from — oo 
to oo. The physical reason behind this 'decompactifi- 
cation' of the angular variable can be understood in the 
following way. In the low energy limit, the gauge field be- 
comes more and more ineffective in scattering fermions 
from one momentum to another momentum along tan- 
gential directions to the Fermi surface. This is because 
the momentum of the gauge field is scaled down under 
the scale transformation, while the circumference of the 
Fermi surface is unchanged. This effectively makes two 
momentum points on the Fermi surface more decoupled 
from each other at lower energies. In other words, the 
'metric' of the Fermi surf ace alo ng the tangential direc- 

e. 




FIG. 5: The shaded areas with angle 6\ and 62 repre- 
sent momentum points which are included in a q (r, x,9\) and 
a q (r, x,92)- The dark region at the center indicate momentum 
points which contribute to both a q (r, x,9i) and a q (r, x, 62)- 
The ratio of the area of the dark region to the area of a strip 
becomes zero in the low energy limit where the momentum 
cut-offs scale as A — > 6~ 2/3 A and A -> &~ 1/3 A. 



along with Eq. (j3"2")h As we go to lower energy (6 > 1), 
the range of the 9 integration increases from (— 7r/2, n/2) 



tions diverges in the low energy limit compared to the 
'metric' along the perpendicular directions. 

Introducing 1+1D fields in real space, 



^ja(r,x,i 



2- 



dw 



A/2 

-A/2 
A/I?l 



dk e l(uT+kx) ij JS (uj,k,9), 



1 /- 00 refill 

a q ( T ,x,9) = — / du d4>\q\ e^ T+x ^a(Lu,q,9- 

2n J-00 J-A/\q\ 



(34) 
(35) 
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we can write down the low energy effective action in the 1+1D real space as 
S = 2_. I drdxdO ip* s (r, x, 9) [d T — isd x ]i(jj S (r, a;, 9) . 

S=±l"' 

+ 2A~ J dTdxded i l«K( r ' x > °) ( 2^ (_a " + q2) + K ) a ^ X > e) 

- . 1/2 ^ / drdxd9dq sa q (r, x, 9)ip* s (r, x,9 + sq/2k F ) ipj s (r, x, 9 

) LI J 

I 



(2tt)i 



sq/2k} 



(36) 



s=±l' 



It should be emphasized that in this action 9 is a non- 
compact variable which runs from — oo to oo. InEq. (|35|) . 
a q (r,x,9) depends on four variables while a(cu,q,9 — (p) 
depends only on three independent variables. The vari- 
able x has been created by the Fourier transformation of 
a(u>, <?, a) with respect to a which is centered at 9. Con- 
ceptually, this is similar to creating a wave-packet which 
is localized in both real space and momentum space by 
linearly superposing wavefunctions whose momenta are 
centered at a particular momentum. The factor 1/2A 
in the second line of Eq. (f36|) is to cancel the double 
counting of momentum points. Note that a q {r, x, 9\) and 
a q {T, x,9 2 ) are not completely independent for different 
9\ and 9 2 . Both a q {T,x,9\) and a 9 (Y, x,9 2 ) include con- 
tributions from a common region in the momentum space 
as is shown in Fig. [5] However, the overlap is not impor- 
tant in determing the scaling dimension of instanton as 
will be shown in the following. The area of the common 
region for a q (r, x, 9\) and a q (r, x, 9 2 ) (the dark parallelo- 
gram in Fig. O is (A 2 /A9) for a small A9 = 9 2 - 9 X . The 
ratio of this area to the area included in a q (r, x, 9) (the 
long strip in Fig. [5]) is 7 = (A 2 /A6»)/(AA). As the mo- 
mentum cut-off decreases as A — > A6~ 2 / 3 , A — ► A6~ 1//3 , 
the ratio decreases as 7 — > 6 _1 / 3 (A/AA0). The ratio 
becomes zero in the low energy limit for any nonzero 
A(9, which implies that the two fields which have a fi- 
nite angle difference (before scaling) arc independent in 
the low energy limit. In the rcscalcd angular variable 
9 = b 1 / i 9, a fixed A9 corresponds to a successively re- 
duced A9 = b~ 1 / i A9 as a low energy limit is taken. 
Since 7 goes as b° for a fixed A9 in the rescaled vari- 
able, two fields which have a fixed angle different A9 
have a finite ratio 7 in the low energy limit. On the 
other hand, two fields whose angle difference increases 
faster than A9 = b° in the rescaled variable have van- 
ishing overlap in the low energy limit. Since the interval 
of 9 increases as b 1 ^ 3 ^, there are infinitely many indepen- 
dent fields a q (r, x, 9) which are separated in the angular 
direction and have only local interactions. We call this 
property an 'asymptotic locality'. This will play a crucial 
role in determining the scaling dimension of instanton as 
will be discussed later. 

Now we can determine the scaling dimension of the 
instanton operator. At a scale set by r and x, fermion 
modes whose angular momenta are less than n max ~ 
fci?min(x, r) = kpx are twisted (at sufficiently large dis- 



tances, we always have r > x because space has the 
smaller absolute scaling dimension). In other words, 
fermion fields that are twisted at scale x have 'wave- 
lengths' larger than 69 ~ l/n max ~ (kpx)" 1 in the space 
of 9. Since x and 9 scale as x — b 2 ^ 3 x and 9 = b~ 1 / 3 9 , 
we have 89' - 6 _1 / 3 (fc f a/) -1 . Note that 59' -> as 
b — > 00 and fermion fields with arbitrarily small 'wave- 
lengths' are twisted in the low energy limit. Since an 
instanton twists fermions of all angle, an instanton cor- 
responds to an operator which creates a vortex with flux 
7r along the non-compact angle direction 9. The physi- 
cal reason why an instanton which is localized in space 
and time becomes an extended object is as follows. In 
the low energy limits, only small angle scatterings are 
important and fermions rarely change the directions of 
their motions. They essentially move on 1+1D planes in 
space and time. Since the distances from the instanton 
and the planes on which fermions move are fixed, at a 
sufficiently large distance scale, all the fermions acquire 
phase 7r as they are transported around the instanton. 
The 1+1D fermions are parametrized by the direction 
of their velocities, 9 and an instanton becomes a vortex 
which is extended along the angular direction. 

It is noted that even though the (Euclidean) Lorentz 
symmetry is broken by the non-relativistic fermions, the 
twist angle will be ir in the long distance limit if there 
is the spatial rotational symmetry as discussed in Sec. 
III. If there is no spatial rotational symmetry, the twist 
angle will depend on 9. But an argument similar to the 
one provided at the end of Sec. Ill can be made to extend 
the conclusion of the following discussion to more general 
cases. In the following, we will focus on the rotationally 
symmetric case. 

The scaling dimension of the extended twist operator 
can be obtained following the reasoning which is analo- 
gous to the state-operator correspondence in relativistic 
quantum field theories. In relativistic cases, instanton 
corresponds to an operator defined at a point in space and 
time. The operator defines a quantum state on the sphere 
enclosing the instanton operator. The scaling dimension 
corresponds to the 'energy' of the quantum state associ- 
ated with the 'time'-evolution in the radial direction. In 
the present non-relativistic case, instanton corresponds 
to an extended vortex operator with flux 7r because in- 
stanton twists boundary conditions of all fermions which 
are parametrized by the non-compact variable 9. The ex- 
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FIG. 6: 'Time '-evolution of a quantum state defined on the 
surface of a pipe extended along the angular direction, where 
a 7r-vortex is pierced through the pipe. Under the time- 
evolution, a point on the surface V is mapped to a point 
on the surface V . 



tended operator defines a quantum state on the surface 
of a pipe V which is extended in the 9 direction in the 
space of r, x and 9 as is shown in Fig. [6] In the func- 
tional Schrodinger picture, the vortex operator defines a 
quantum state as 



m> js (T, x, 9),<ip* s (T, x, 9), a q {T, x, 9)] 
= [ DiP J8 DiP* a Da q e- s ^*' a l 



(t,x,6)&V 



(37) 



On the r.h.s. of the above equation, the fermion fields 
have the anti-periodic boundary condition and all the 
fields inside the pipe are integrated out with the condition 
that the fields on the surface of the pipe coincides with 
the fields tpj S {r, x, 9) 7 ip* s (T, x, 9), a q (r, x, 9). The scaling 
dimension corresponds to the 'energy' of this quantum 
state associated with the time evolution given by 



x 

T 



b 2 / 3 x, 

6V, 

b-^ 3 9. 



(38) 



We can define a 'Hamiltonian' for the 'time' evolution 
because the action is local in r and x. Namely, a quantum 
state on a surface with larger \t | and \x | are uniquely 
determined from the state on a surface with smaller \t\ 
and \x\. 

Since r and x scale differently, the surface of the pipe 
should 'expand' in different rates depending on its normal 
vector. Note that the 'time' evolution also involves the 
transformation in 9 because of the anomalous dimension 
of the angular variable. In determining the scaling di- 
mension of instanton, the key is the locality of the action 
in the angular variable 9. The action Eq. (|36p is asymp- 
totically local in the space of 9 in the following senses. 



First, angles of the fermions can change at most by A/fcp 
through an interaction with the gauge field. Second, the 
fermion field at angle 9 is coupled only with the gauge 
field near angle 9. Third, two gauge fields which have dif- 
ferent angles 9\ and #2 are independent for a sufficiently 
large \9i — 9\\ which is, yet, much smaller than the range 
of the 9. One may worry about a possible breakdown of 
the locality in the angular direction in the presence of 
short range four fermion interactions. Indeed, the four 
fermion interaction 

V I Uf^ckJidkidei tp* (wi, fci, W(w 3 , h, 9 3 ) 
i>(u2, hi, 2 )V(w 4 , hi, 9 i )6(J2(-iy^) 

i 

*($3(-1)*(*f + h) cos 0i)6(52{-lY{k F + hi) sin ft) 

i i 

(39) 

is non-local in the angular direction. However, the inter- 
action strength scales as V — > b^ 2 ^V under the trans- 
formations in Eqs. (J32)) and (|33[) and the four fermion 
interaction is irrelevant at low energies. 

Because of this locality of the action, any extended 
object should have an 'energy' which is either zero or 
infinite with respect to the vacuum. This is a very gen- 
eral statement for a local theory. For example, if a local 
theory is defined in a range < 9 < 2L where L is 
much larger than any length scale of the local coupling 
in the theory, the energy of the system is roughly the 
twice of the energy of the system defined in the range 
< 9 < L. This implies that only or ±oo are pos- 
sible for the energy of the vortex operator with respect 
to the vacuum energy. The vortex is a non-trivial object 
which necessarily 'excites' the fcrmionic state by twist- 
ing the boundary condition and it should have an infinite 
'energy'. Note that — oo (0) are excluded because a neg- 
ative (zero) energy would imply that instanton becomes 
more relevant (equally relevant) with increasing number 
of spinon flavors or with increasing length of Fermi sur- 
face in the momentum space. This is unlikely because 
the fermions always screen the gauge field and the scal- 
ing dimension of instanton should increase as the number 
of fermion modes increases. Although the scaling dimen- 
sion defined in the thermodynamic limit is infinite, the 
scaling dimension of the regularized instanton operator 
is well defined and systematically increases as the num- 
ber of available fermion modes increases. The number of 
fermion modes in a finite system is proportional to the 
number of spinon flavors and the length of the Fermi sur- 
face due to a finite mesh in the momentum space. This 
implies that the scaling dimension of the instanton oper- 
ator diverges more rapidly as the long distance limit is 
taken when there are more flavors or longer Fermi sur- 
face. An explicit example of this is provided in Eqs. (fT6|) 
and (|17p for the non-interacting system, where the scaling 
dimension of the regularized instanton operator increases 
as N or kp increases. Therefore we conclude that instan- 



ton has to have a positive infinite scaling dimension 
general. 



in 



general 



V. STRONGLY INTERACTING FERMIONS 
WITH A FEW FLAVORS 

Now we move on to the physical case with a small 
but nonzero N, e.g., N = 2. In this case, even if 
one can somehow ignore instantons, the fermions are al- 
ready strongly coupled with the non-compact gauge field. 
Therefore, one can not exclude the possibility that the 
strongly interacting fixed point becomes unstable against 
a particle-hole or particle-particle condensation 21, 37, 
[38[ due to the strongly fluctuating non-compact gauge 
field. Here we set aside those possibilities caused by 
the non-compact gauge field and focus on the question 
whether the fixed point is stable against proliferation of 
instantons or not. Because we can not ignore vertex cor- 
rections any more, we do not know what the precise form 
of the scale transformation is for the strongly interacting 
fixed point. However, we can still argue that the scaling 
dimension of an instanton is infinite based on the locality 
of the low energy theory in the angular direction. The 
emergence of the locality in the angular direction is in- 
dependent of a particular form of scale transformation. 
Because of the Fermi surface geometry, momentum of the 
gauge field should be scaled down more slowly compared 
to radial momentum of the fermions. If a momentum 
of fermion scales as k — b~ a k then a momentum of the 
gauge field should scale as q = b~ a / 2 q , which is the con- 
sequence of the fact that Fermi surface is locally parabolic 
unless there is a nesting which we do not consider here. 
Therefore, we have q/k >> 1 at low energies and this 
forces fermions at a certain angle arc coupled only with 
the gauge field whose momentum is tangent to the Fermi 
surface. This guarantees that the low energy effective ac- 
tion should be local in the angular direction. Moreover, 
the mixing between k and q gives rise to the anomalous 
scaling dimension for the angular variable, a/2 > and 
the angular variable becomes decompactified in the low 
energy limit. Again, this implies that the instanton oper- 
ator is an extended vortex with flux tt along the extended 
angular direction. Because of the locality, the extended 
vortex should have an infinite scaling dimension as we 
have discussed previously. 



VI. CONCLUSION 

In conclusion, we argue that the U(l) spin liquid 
state is stable against proliferation of instantons for any 
nonzero number of spinon flavors if there is a spinon 
Fermi surface. We formulated the low energy modes 
near the Fermi surface in terms of an infinite set of 
chiral fermions and made an observation that the an- 
gular variable that parametrize the Fermi surface ac- 
quires a positive scaling dimension due to quantum ef- 
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fects and becomes a non-compact variable in the low en- 
ergy limit. Because the low energy effective theory is lo- 
cal in the non-compact angular direction, an instanton, 
which twists boundary conditions of all chiral fermions, 
should have an infinite scaling dimension. Therefore, 
instantons are strongly irrelevant and the non-compact 
U(l) gauge theory is a good low energy description. 

A few comments are in order. First, since the scaling 
dimension of instanton is infinite in the infrared limit, a 
small fugacity of instantons will rapidly flow to zero at 
low energies. A finite fugacity at an intermediate length 
scale implies that there is a small but finite density of 
instantons at the length scale. However, the renormal- 
ization group flow implies that the density at a larger 
distance scale will become smaller and instantons are not 
important in the long distance limit. If the fugacity y is 
tuned to a sufficiently large value via tuning some micro- 
scopic parameters, then the nonlinear terms in the flow 
equation, Eq. ()18|) becomes important and the sum of 
them may not converge. This signifies a breakdown of 
the perturbative expansion which is valid within a finite 
domain near the y = point. If this happens, the fu- 
gacity of the instanton operator can flow toward a large 
value leading to a confinement, despite the fact that the 
linear term has a large negative coefficient. It is noted 
that the infinite scaling dimension of instanton at the de- 
confined fixed point y = does not necessarily imply that 
the confinement phase is always unstable. Both the de- 
confinement phase and the confinement phase may have 
finite regions of stability in the parameter space, sepa- 
rated by a phase transition. To understand the nature of 
the confinement phase is an important open problem [4fJ. 
Second, we expect that the present argument can be gen- 
eralized to the cases where there is no rotational symme- 
try or there arc only segments of Fermi surface. This 
has been already demonstrated for the free fixed point at 
the end of Sec. III. We expect that the similar argument 
will hold true at the interacting fixed point too. This is 
because any finite segment of Fermi surface contains an 
infinite number of modes which contribute to the scaling 
dimension of instanton. 



VII. ACKNOWLEDGMENT 

This work has been supported by NSERC. The author 
thanks Matthew Fisher, Yong Back Kim, Subir Sachdcv 
and Xiao-Gang Wen for helpful discussions and, partic- 
ularly, Patrick Lee for helpful comments and suggestions 
to improve the paper. 



APPENDIX A: Scaling dimension of a general twist 
operator 

Here we calculate the scaling dimension of a general- 
ized twist operator. Suppose the boundary condition of a 
fermion is twisted by an arbitrary angle around the twist 
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FIG. 7: Feynman diagrams for the one-loop self energies of 
(a) the gauge field and (b) the fermions. The solid (wiggled) 
line represents the propagator of the fermions (gauge field). 



operator which is inserted at the origin. The fermion 
satisfies the twisted boundary condition 

^(ze 2m ) = e J< V(z) (Al) 

in the complex plane with z = r — ix, where $ = 
27r(l/2 — 77) with < 77 < 1. 77 = 1/2 corresponds to 
the untwisted periodic boundary condition and r\ = 0, 
the antiperiodic boundary condition. A general value of 
r\ between and 1 corresponds to a fermion twisted by 
a phase angle between — 7r and 7r. The twisted boundary 
condition is implemented by the mode expansion, 



<l>(z) 



E 







n+W 



z n+ri+l/2 '■ 



2-^1 7 n~r 1 +l/2 1 



(A2) 



where i/j n+r > and "0*™ -r ? with integer n are the n-th nor- 
mal modes. Note that ip an d ij)* have the different mode 
expansions because the gauge invariant operator %j)*ip al- 
ways have to satisfy the untwisted periodic boundary 
condition. The scaling dimension of the twist operator is 
given by 



dz , , 

17- zT(z) 



(A3) 



twisted b.c. 



where the contour of the integration encloses the origin 
and 



T(z) = I : [(cty*)V - V*dV] : 



(A4) 



is the regularized holomorphic energy-momentum tensor 
with 



: if)* (z)dip(z) := lim 



i/j*(w)dip(z) 



1 



(w — z)' 



(A5) 



To calculate < T(z) > twisted b.c, we first consider the 
expectation value of a fermion bilinear with the twisted 
boundary condition, 



( V>*HV>(z)> twisted b.c. 



= E 



yjm+l/2-rj z n+l/2+ri 
Z 



n .-n-\/2 -n+1/2 



w 



(A6) 



where we have used the commutator {ip m ^ ,ip n+rl } = 
Sm t - n and the property of the vacuum, ip n+v \0 >= 
^™-'?t|0 >= for n > OHI. From Eqs. |A4l) - (lA6J) . 
we obtain 



*C 1 \ZJ ^twisted b.c.~ 



(27? - 1) 



(A7) 



and from Eq. (|A3|) the scaling dimension is obtained to 

be 7 • In the presence of multiple fermions, the total 
energy momentum tensor is the sum of individual energy 
momentum tensors and the scaling dimension is the sum 
of individual scaling dimensions. If the i-th fermion is 
twisted by an angle $ j , the scaling dimension of the twist 
operator becomes 



^E 



$? 



r2' 



(A8) 



APPENDIX B: Calculation of the self energies 

In the large N limit, the vertex corrections are negli- 
gible and the one-loop corrections are dominant [191]. Al- 
though the one-loop self e nerg ies of the fermions and the 
gauge field are well known [19l |20|, l2l| . it is instructive to 
calculate them from the effective action Eqs. (|27[) - ([2"5)) 
to make sure that the theory contains the essential low 
energy physics. 

From the quadratic action in Eq. (|27p , the propagator 
of the fermion is given by 



^.(w.MWvK*,* ))=<*,■,*'*. 



3 s 



3,3 s,s 



5(uj - lo )8{k - k )6(0 - 9 )g s (u>, k, 0), (Bl) 

where g s (ui,k,8) = iuJ } sk , and the propagator of the 
gauge field is given by 

a(u, q, a)a{v ,q ,a )) = y- t5{v + v )S(q + q )S(a — a ) 



T> (u!,q, a), 



(B2) 



where T>q{v, q, a) = 1 . 2 } 2\<n K - The unconventional 

factor, 1/|<7| in the gauge propagator is due to the angular 
representation. 
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Applying the standard Feynman rule to the action in 
Eq. ([25]). we can calculate the self energy of the gauge 
field (Fig. [3(a)) as 



IL(v,q,a) 



J^Yl J <^dkdcf> 



g s (uj + v/2, k + q<f)/2, (f> + a + sq/2k F ) 
g s (uj — v/2, k — q<f>/2, <fi + a — sq/2kp) 
Nk F Nk F \v\ 



2ir A-k \q\ 

for |g| >> \v\. The dressed gauge propagator becomes 



(B3) 



V-\v,q, 



2K 



Nk F 
~2^ 



7l 



Vi\ 



Xq\ (B4) 



where 7 and \ are constants. The first term should be 
canceled due to the gauge invariance. The diamagnctic 
term K depends on the details of the high energy cut-off 
(rcgularization) scheme and we fix its value by requiring 



the gauge invariance. From Fig. [7](b), we obtain the self 
energy of the fermions as 

£ s (w, k, 6) = -t^t^ / dvdq\q\d(j> 

T>{v, q,9 — (p)g s (u> — v, k — q<fi, — sq/kp), (B5) 

where it is essential to use the dressed gauge propagator, 

1 



V{y, q, a) 



7IH/M +XQ 2 



(B6) 



The integration can be done straightforwardly and we 
obtain, 



£ 5 .( W ,fc,0)= iC sgnMM 2 / 3 , 



(B7) 



where c ~ -y -1 / 3 ^- -2 / 3 is a constant. Even though we use 
the dressed fcrmion propagator, the leading behaviors of 
the self energies are not modified [19|. 
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